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Abstract 

The Debye screening masses of the a, lo and neutral p mesons and the photon are calculated in 
the relativistic mean-field approximation. As the density of the nucleon increases, all the screening 
masses of mesons increase. It shows a different result with Brown-Rho scaling, which implies 
a reduction in the mass of all the mesons in the nuclear matter except the pion. Replacing the 
masses of the mesons with their corresponding screening masses in Walecka-1 model, five saturation 
properties of the nuclear matter are fixed reasonably, and then a density-dependent relativistic 
mean-field model is proposed without introducing the non-linear self-coupling terms of mesons. 

PACS numbers: 21.65.+f, 13.75.Cs 
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I. INTRODUCTION 



In spite of the great success of the relativistic quantum field theory(RQFT) in many 
areas of modern physics, some difficulties have been found in applying this theory to nuclear 
systems. In nuclear physics, the coupling constant of strong interaction between nucleons is 
far larger than the fine structure constant in quantum electrodynamics(QED). The ground 
state of the nuclear matter or the finite nuclei is often defined as "vacuum" , where the Fermi 
sea is filled with interacting nucleons, and no anti-nucleons and "holes" exist. 

In 1970's, in order to solve the nuclear many-body problems with RQFT, Walecka et 
al. developed the quantum hadrodynamics theory (QIID)[0, ^. It is a important progress 
in nuclear physics. Since then, the QHD theory has widely been used in nuclear physics 
and nuclear astrophysics [0, However, in the earliest QHD theory, the resultant com- 
pression modulus is almost 550MeV|l[], which is far from the experimental data range of 
200 — SOOMeV. To solve this problem, the nonlinear self-coupling terms among a mesons 
were introduced in addition to the mass term \m'^cr'^ 0, 

U{a) = \mla' + \g,a' + \g,a\ (1) 

Moreover, a self-coupling term of the vector meson 

is added to produce the proper equation of state of nuclear matter]^, |^. No doubt, addi- 
tional parameters would give more freedoms to fit the saturation curve of nuclear matter. 
Zimanyi and Moszkowki developed the derivative scalar coupling model yielding a compres- 
sion modulus of 225MeV without any additional parameter P]. Because all of these models 
include the nonlinear density-dependent terms in the Lagrangian in substance, the proper 
value of the compression modulus can be obtained. 

In the relativistic quantum field theory, the nonlinear density-dependent terms in the 
Lagrangian represent higher-order corrections. In this paper, we would calculate the self- 
energies of virtual mesons with zero time-space momentum in the framework of the rela- 
tivistic mean-field approximation, and then the Debye screening masses of the mesons would 
be discussed. Finally, we would replace the masses of the mesons with their correspond- 
ing screening masses in the Walecka-1 model, and develop a density-dependent relativistic 



2 



mean-field theory without additional parameters, which is similar to the density-dependent 
quark-meson coupling model, i.e. quark-meson coupling model-2(QMC-2)[0. 



II. THE SCREENING MESON MASSES IN THE RELATIVISTIC MEAN-FIELD 
APPROXIMATION 

the Lagrangian density in nuclear system can be written as 

£ = Cn + Ca + + Cp + Lint- (2) 

In this expression, the Lagrangian density for the free nucleon field can be described by: 

Ln = ^{i-f,d^-M^)^, (3) 

where ip is the field of the nucleon and M^r is the bare mass of the nucleon. The free 
Lagrangian densities for the a, u, p meson fields and the photon field can be expressed by: 

= \d,ad^a - \mla\ (4) 

l^u. = —uJ^uUJ^'' + ^mluj^uj^, (5) 

= -\r,.-R^'' + \mlp,-f?^, (6) 

C-A = -\f,^F^\ (7) 

where 

i-^f,u = df,ujy - duUJ^, (8) 

Ri,u = d^p^ - dyp^,, (9) 

= d^A,-d,A^ (10) 

are corresponding field tensors, respectively. The interactive Lagrangian density can be 
written as 

- gp^l'l ■ P> - eij^.^^A'^^lj (11) 

with r being the Pauli matrix. 
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In the framework of relativistic mean-field approximation, the meson fields operators and 
the electromagnetic field operator can be replaced by their expectation values in the nuclear 
matter 0: 

(12) 
(13) 
(14) 
(15) 

where uo is the expectation value of the scalar meson field operator, c^o and po are the 
expectation values of the time-like components of uo meson and neutral p meson field oper- 
ators, respectively, and Aq denotes the scalar potential of the electromagnetic field in the 
nuclear system. The Lagrangian density for nuclear matter in the relativistic mean-field 
approximation then reads 

^RMF = i' iil^id^ - Mn) tp - Qa'ipipaQ - g^ip-f^ipujo - gpipj^—ippo 

^ 22,-'- 22,-'- 22 

In the nuclear matter, because of the long-range Coulomb repulsive interaction between 
protons, the nuclear matter would not be bound. As a model of nuclear matter, the electro- 
magnetic interaction usually be ignored. However, in the finite nuclei, the electromagnetic 
interaction should be considered in the structure calculation. 

The screening mass of the meson m* is defined from the static infrared limit = and 
then — i> as Ref. ||10 

^Jml + S„(/c = 0,pp,p„), (a = a,u;,po). (17) 



It represents the inverse Debye screening length and implies the long-distance correlations. 
In the following, we will calculate the self-energy of virtual mesons with A; = in the 
framework of relativistic mean-field approximation, and then the screening masses of the 
mesons and the photon are obtained. 

The interactive perturbation Hamiltonian of the model can be expressed as 

Hi = g^ija^j + g^^J-ipU^^^J + gpi^i^^ ■ Pp^ + e^j^iJ^^^A^^, (18) 
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and the S-matrix can be written as 



S ^ So + Si + S2 + (19) 

where 

Sn = j d^xi j d^X2 ... j d^XnT [Hi{xi)ni{x2) . . . (20) 

In the 2nd-order approximation, only 

-^2 = ^ j d'x, J d'x2T[ni{xi)ni{x2)]. (21) 

should be calculated. 

The nucleon field operator ip{x) and its conjugate operator ■0(x) can be expanded in 
terms of a complete set of solutions of the Dirac equation: 



[ApxUr,{p, A) exp [ip-x- ie^'^\p)t) 
+ Bl^Vr,{p,X)exp{-zp-x-ze^-\p)t)], (22) 



|^Aj^;^f7,,(;), A) exp {-ip- X + ie^~^\p)t) 
+ BpxVr,{p, A) exp {ip-x + ie^-\p)t)] , (23) 

where Urj{p, A) and K?(p, A) are Dirac spinors for the positive and negative energies, respec- 
tively, and 

J][/,(p,A)[7,,(p,A) = (24) 

A=l,2 ^ ^ ^ 

^K,(p,A)K,Kp,A) = (25) 

A=l,2 ^ ^ ^ 

with A and rj denoting the spin and isospin for the nucleon, and — + g^ao is the 
effective mass of the nucleon. 

£(±)(p) = ±E*{p) + g^uo + gp^po + e^^Ao (26) 
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with E*{p) = yjp + M]^^, are the positive and negative energy eigenvalues for the Dirac 
equation of the nucleon in the nuclear matter, respectively. Assuming there are no antinu- 
cleons in the nuclear matter or finite nuclei, only positive-energy components are considered 
in Eqs. (|2|) and (||). 

When a scalar meson of momentum k is considered in nuclear matter, the scalar meson 
field operator a{k,x) can be expressed as 

a{k^x) = a{k) exp{—ik ■ x) + a\k) exp{ik ■ x). (27) 

If the coupling constants g^, Qui and Qp, the proton charge e, and the masses of mesons 
are supposed have already been renormalized, the contribution of a single nucleon loop is 
not necessarily calculated [|Il|], only these contributions from the Feynman diagrams in Fig.l 
should be considered. 

The expectation value of 5*2 can be written as 

{k2\S2\k,) = - igl{2'K f5\pr + k, - - k^) 



r;=l,2 A=l,2 



where ki = k2 = k, and Pi = P2 = P = {E*{p),p), and 9{x) is the step function. 

Considering the diagrams in Fig.l, the scalar meson propagator G{k) in nuclear matter 
can be derived as 

1 i 1 i 



(27r)4 F - m2 + is (27r)4 p - + is 

?7=1,2A=1,2 -J \ J yfJ 



i)-^-M*j^ i) + ^-M*^J ' (27r)4 P - m2 + i£ 

I 1 I gl ^ ^ f d^p M*^ 



(27r)4 p - m2 + (27r)4 - m2 + _ 
According to the Dyson equation 



we obtain the self-energy of the scalar meson in the nuclear matter 

d^p Mil 



9^ E E 

J7=l,2 A=l,2 



£/,(p,A) 



+ 



N 



(31) 



where 



A — 1)2 



1 



M 



A=l,2 

Tr 



N 



+ 



1 



f/,(p,A) 



N 



U,ip,X)U,ip,X) 



{p-ky-M*^^ ' {p + k)^-M*^ 



2M^ J 



(32) 



In the above deduction, Eq. is used. 

In order to obtain the Debye screening mass of the scalar meson, we calculate the self- 
energy of the scalar meson in the nuclear matter in the static infrared limit = and then 
k 0. Because the time-space momentum of virtual scalar meson is zero, /c = 0, we have 
k^ = 0. Considering = Mj^^, the virtual scalar meson self-energy 

d^p M*^ 



E. 



9l E 

r?=l,2 



(27r)3E*(j9) 



0{PF-\P\] 



Tr 



i>-f: + M'„ i> + f: + Ml,\'l,+ M-^ 



A E 

»7=1,2 

^% (2vr) 



—2p ■ k 2p ■ k 

d^p Mil 



2M 



N 



(27r)3E*(p) 



(^{Pf - \p\)Tr 



N 



p ■ k 2M* 



N J 



M* 



(33) 



where and are the scalar densities of protons and neutrons with 

d^p 



p 
Ps 



Ps 



Mm 



d^p 



Mil 



n {2rry ifP + Mli^)-2 



(34) 
(35) 
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Similarly, the self-energies of the to", neutral p and 7 with /c = may be obtained: 

^ _ 91{PI + P's) 



2M^ 



(36) 



" 8M^ ' ^^^^ 

2 p 

= li^- ^^^^ 
According to Eq. (0), the Debye screening masses of the a, uj neutral p and the photon 
can be written as 



m 



ml{pN,fp) = ^/ml + S^, (39) 
:{PnJp) = V< + (40) 



m*n{pNjp) = \ ml + Sp, (41) 



^a(p^,/p) = V^^ = \Ij^- (42) 

They are also the functions of the density of nuclear matter p^r and the proton fraction 
fp = Z/A. 

In the relativistic mean-field approximation, 

mlao = -gM + Ps)^ (43) 
so the effective mass of the nucleon can be written as 

M*^ = Mn + ga(Jo 

= M, - 3lif^±M. (44) 

From Eq. (^), the screening mass of the scalar meson can be expanded as 

1 9l P's +_Ps 



To the first order of the scalar density of the nucleon, 

where 6mcr = m* — rricr, and 6Mn = — Mn, At the low density of the nuclear matter, 
Ml^K. Mi^, and 

5m^ 5M^ _ _1 
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Obviously, the fractional change of the mass of the scalar meson is -1/2 times that of the 
nucleon in the nuclear matter, and this simple value of the ratio is similar to the scaling 
relation in Ref. 

The ratio of the effective mass to the bare mass for the nucleon and the ratios of screening 
mass to the bare mass for mesons as the functions of the nucleon density for the symmetric 
nuclear matter are shown in Fig. 2. It is seen that when the number density of nucleon 
increases, the effective mass of the nucleon decreases, and the screening masses of mesons 
increase. The increase of screening masses of the mesons means that as the density of the 
nuclear matter become dense, the range of nuclear forces decreases, and the Debye screening 
effect enhances in the denser nuclear matter. Our results on the screening masses of the 



mesons are different from the Brown- Rho scaling[|l^, in which the masses of the mesons are 
reduced in the nuclear matter, similarly to the reduction of the effective mass of the nucleon, 
as the density of the nuclear matter increases. 



III. THE EQUATION OF STATE FOR NUCLEAR MATTER 

The relativistic mean-field results may be derived by summing the tadpole diagrams self- 
consistently in nuclear matter, retaining only the contributions from nucleons in the filled 
Fermi sea in the evaluation of the self-energy and energy density]^], so the relativistic mean- 
field approximation is consistent to the relativistic Hartree approximation in the calculation 
of nuclear matter. Because the time-space momentum k of virtual mesons in the tadpole 
diagrams is zero approximately, if the Debye screening effect is considered, the masses of 
mesons in the relativistic mean-field approximation should be replaced by their screening 
masses, respectively. Therefore, the total energy density and the pressure of nuclear matter 
can be deduced to 

1 *2 2 1 *2 2 -'- *2 2 , i'AO\ 

^=2^ -^0 - 2^ co^o - 2^ pPo + 2^ ^B, (48) 

B=p,n 

and 

B=p,n 

1 *2 2 I *2 2 I -'- *2 2 f An\ 

-Tj^ <x^0 + 7^^ u^Q + pPo^ (49) 
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with 



es = —-3 dp[{f + M*^y + (g^u, + gp^Po)) , (50) 



(2vr) 

and the vector densities of protons and neutrons being 



respectively. 

By fitting the saturation properties of nuclear matter, the parameters of the relativistic 
mean-field approximation in which the mesons have density-dependent screening masses can 
be fixed 



0^ g2 

^ = 8.297/m^ ^ = 3.683/m^ ^ = 5.187/ml (53) 

rri'^ 



With these parameters we obtain a saturation density of 0.149fm~^, a binding energy 
of 16.669 MeV, a compression modulus of 280.1 MeV, a symmetry energy coefficient of 
32.8 MeV and an effective nucleon mass of O.8O8M7V for the symmetric nuclear matter. 

The saturation curves for the nuclear matter with different parameter sets are plotted 
in Fig. 3. Comparing these results with each other, we see that when the density of the 
nucleon becomes larger than the saturation density of nuclear matter, the average energy 



per nucleon with our model parameter set increases more slowly than those with NL3[|13 
NLSH[0 and TM1|]^ . It manifests that the equation of state for nuclear matter in our 
model is softer than those in the other models, because the mesons have density-dependent 
screening masses. 

Due to the approximate chiral symmetry restoration in nuclear matter, the effective 
nucleon mass is reduced. Brown-Rho scaling would imply a similar reduction in the mass of 
all the mesons except the pion. However, this scaling would not lead to reasonable equation 
of state for the nuclear matter. Our model gives different results with Brown-Rho scaling, 
and can fit the correct saturation properties of nuclear matter in the framework of relativistic 
mean-field approximation. 

In the nuclear matter. The screening mass of photon is only related to the scalar density 
of the proton, but not the momentum of the photon. In a symmetric nuclear matter, if 
one takes a nucleon number density of 0.16/m~^, and the bare nucleon mass, the resultant 
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screening mass of photon is about 5A2MeV [0. Although the photon gains an screening 
mass in the nuclear matter, the range of Coulomb repulsive force is still large enough, and 
the nuclear matter would not be bound as the electromagnetic interaction is considered. In 
the case of finite nuclei, the contribution from the photon mass term should be included in 
the relativistic Hartree approximation or relativistic mean-field approximation. 

In nuclear matter, the screening masses of the mesons increase with the density of the 
nuclear matter. This result is equivalent to the statement that the coupling constants 
decrease with the increasing number density of nucleon when the masses of mesons retain 



constant. At this point, our model is consistent to the model in Ref. [16 



IV. SUMMARY 



In summary, we calculate the self-energies of virtual mesons with zero time-space mo- 
mentum in the relativistic mean-field approximation, and the Debye screening masses of 
the mesons are obtained, which are the functions of the scalar densities of protons and neu- 
trons. The screening masses of the mesons increase with the density of the nucleon. It shows 
different results with Brown-Rho scaling. Replacing the masses of mesons with correspond- 
ing screening masses in the relativistic mean-field approximation, we obtain a relativistic 
density-dependent nuclear model. In this model, the nonlinear self-coupling terms of mesons 
are not needed and only three model parameters are required. With this model, five satu- 
ration properties of symmetric nuclear matter, the saturation density, the binding energy, 
the effective nucleon mass, the compression modulus and the asymmetric energy, are calcu- 
lated. In denser nuclear matter, the equation of state with this model is softer than those 
of previous models. This implies that the dense matter in the core of neutron stars or the 
center of the relativistic heavy ion collision might be described correctly. 
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KJCX2-N11 and the National Natural Science Foundation of China under grant numbers 
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1 3 



FIG. 1: Feynman diagrams for the meson or the photon self-energy in nuclear matter, while 1 and 
2 denote particles of the initial state, 3 and 4 denote particles of the final state. 
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FIG. 2: The ratios of the effective or screening mass to the bare mass M* /M for the nucleon and 
the mesons as a function of the nucleon density /J at in the symmetric nuclear matter in this model. 




FIG. 3: Average energy per nucleon e/ pN — Mn as a function of nucleon density pjy with the 
different parameter sets. The solid line denotes the curve obtained in this model, the dash-dot line 
is for the NL3 parameters |jl3|, the dash line is for the NLSH parameters , and the dot line is for 
the TMl parameters 0. 



14 



